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Abstract 

Intertrade duration of equities is an important financial measure characterizing the trading activities, which is defined 
as the waiting time between successive trades of an equity. Using the ultrahigh-frequency data of a liquid Chinese 
stock and its associated warrant, we perform a comparative investigation of the statistical properties of their intertrade 
duration time series. The distributions of the two equities can be better described by the shifted power-law form than 
the Weibull and their scaled distributions do not collapse onto a single curve. Although the intertrade durations of 
the two equities have very different magnitude, their intraday patterns exhibit very similar shapes. Both detrended 
fluctuation analysis (DFA) and detrending moving average analysis (DMA) show that the 1-min intertrade duration 
time series of the two equities are strongly correlated. In addition, both multifractal detrended fluctuation analysis 
(MFD FA) and multifractal detrending moving average analysis (MFDMA) unveil that the 1-min intertrade durations 
possess multifractal nature. However, the difference between the two singularity spectra of the two equities obtained 
from the MFDMA is much smaller than that from the MFDFA. 
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1. Introduction 



The intertrade duration, defined as the waiting time between two consecutive transactions of an equity, is at the 
core of the autoregressive conditional duration (ACD) model Jll-Ql, which is deemed as the transaction -level version 
of the Nobel Prize winning autoregressive conditional heteroskedasticity (ARCH) model 01. With the availability of 
ultrahigh-frequency data, the statistical properties of the intertrade durations of some stocks have been investigated 
and some common and idiosyncratic patterns have been unveiled. 

The most studied property is the distribution of intertrade durations of different stocks in different markets. There 
is conclusive evidence showing that the survival function of durations cannot be represented by a single exponential 
I5-I8J. Empirical analyses proposed diverse functional forms for the distributions, such as the Mittag-Leffler function 
I9, fl0ll . power laws ifloL 11], modified power laws 1 12, 13], stretched exponentials (or Weibulls) 1 14-3, stretched 



exponentials followed power laws 1 2(| 21 ] . Recently, statistical tests were performed to illustrate that the waiting time 
distributions are neither exponentials nor power laws and the ^-exponentials and the Weibulls are better candidates 
1 2211 . For some stock markets, it has been found that the intertrade durations of different stocks exhibit a scaling 
behavior 1 3, IH and the distribution is Weibull followed by a power law tail lfl9ll . 

The long-term correlation or persistence in intertrade durations plays an important role in the ACD-type models 
1 230 ■ In the econophysics community, the first research was conducted on 30 stocks listed on the New York Stock 
Exchange (NYSE) from January 1993 to December 1996 using the detrended fluctuation analysis (DFA) approach 
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In addition, crossover behaviors have been unveiled in the fluctuation scaling for many stocks on the NYSE, 
NASDAQ and Shenzhen Stock Exchange in different time periods J24l - l26ll . Furthermore, the multifractal nature in 
the intertrade durations of German DAX stocks (from 28 November 1997 to 31 December 1999) and SZSE stocks 
(from 2 January 2003 to 3 1 December 2003) has been investigated and confirmed based on the multifractal detrended 



fluctuation analysis (MFD FA) approach II261I27H 



The brief but complete review presented above shows that the statistical properties of intertrade durations have 
been studied only for stocks and relatively fewer studies have been conducted on the linear and nonlinear long-term 
correlations. In this work, we perform a comparative analysis of the statistical properties of the intertrade durations 
of a very liquid Chinese stock and its associated warrant. Special emphasis has been put on the correlations of the 
durations. Comparing with the data in Ref. |ljl 26|, the stock investigated in this work is listed on the Shanghai 
Stock Exchange (SHSE) and the data were recorded during a very bullish market phase. The results show that the two 
equities share very similar statistical properties. 



2. Data description 

In this work, we analyze the trading records of the Bao Steel stock and its warrant (a kind of call options) traded 
on the SHSE. The records started on 22 August 2005 and ended on 23 August 2006, containing 243 trading days. 
There are more than 25 million trades for the stock and about 3.8 million transactions for the warrant. Each data set 
records the relevant information of all the transactions. The time stamps are accurate to 1 second. For each equity, 
the time series of intertrade durations, denoted r(0, is determined. In this procedure, only trades in the continuous 
double auction (from 9:30 a.m. to 11:30 a.m. and from 13:00 p.m. to 15:00 p.m.) are considered and the waiting 
times between the last trade in the morning and the first trade in the afternoon are excluded from the analysis. We find 
that 91.33% intertrade durations of the stock are null and 82.44% durations of the warrant are zero. 



3. Basic statistical properties 



3.1. Probability distribution 

In order to compare the distributions of the two equities, we investigate the normalized durations t/o~ t , where cr T 
is the sample standard deviation of t lfl9ll . Figure 02a) plots cr T p(T) as a function of r/cr T for the Bao Steel stock and 
its warrant. The most significant feature is that the scaled probability distributions of the stock and its warrant do not 
collapse onto a single curve, which is different from the scaling pattern for stocks IU6U19I1 . 
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Figure 1: (Color online) Empirical distributions of the scaled intertrade durations of the Bao Steel stock and its warrant in the whole time period 
(a) and in four successive time periods (b). 
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Following Refs. 22], we adopt the Weibull and the shifted power-law distributions to model the scaled inter- 



trade durations, g = t/ct t . The Weibull probability density p w (g) can be written as 

Pw {g) = uPff- 1 exp(-a^). (1) 

When p = 1, p w {g) recovers the exponential distribution. When < /3 < 1, p w (g) is a stretched exponential or 
sub-exponential. When /3 > 1, p w (g) is a super-exponential. We adopt the maximum likelihood estimator for model 
calibration and find that a - 1.22 and /3 = 1.41 for the Bao Steel stock and a = 0.75 and /? = 0.97 for the warrant. It 
means that the intertrade durations of the stock are sub-exponentially distributed, while the warrant has an exponential 
distribution in the intertrade durations. The resulting fits are illustrated in Fig. [JJa) as dashed lines. It indicates that 
the Weibull distribution is not suitable for median and large durations. 



The shifted power-law distribution is actually also known as the ^-exponential distribution 112 8l 12911 . which is 
defined by 

Pci (g) = a(g + g Q r^ 1 \ (2) 

where go > 0. When g » go, we observe a power-law behavior in the tail with a tail exponent of y. Nonlinear 
least-squares regressions give that go = 4.09 and y = 5.50 for the stock and go = 3.76 and y = 5.70 for the warrant. 
The resultant fits are illustrated in Fig. [TJa) as continuous lines. We find that the tail behaviors of the durations are 
quite similar. 

The difference in the intertrade duration distributions between the current stock and the SZSE stocks studied 
in Ref. Jl9ll stems from the different time resolutions of the data. The super-exponential distribution describes the 
durations less than one second, while the sub-exponential distribution here concentrates on durations longer than one 
second. For median and large intertrade durations, the distributions are quite similar. 

In order to investigate the stability of the distributions, we partition each data set into four segments with almost 
the same number of trading days. The empirical distributions in each time intervals for the two equities are plotted 
in Fig. [Tib). It is obvious that the four distribution curves for each equity collapse excellently onto a single curve. It 
shows that the distribution is stable during different time intervals. 

3.2. Intraday pattern 

The intertrade durations have been found to exhibit an inverse {/-shaped intraday pattern in the NYSE market 



II lL I3QI1 and Russian stock market Bill , which indicates higher trading activity in the open and close than that in other 
time during each trading day. For the SZSE stocks traded in 2003, a different and more complex intraday pattern has 
been uncovered 12^1 . It is interesting to check if the intraday pattern changed when new data are used. 

We segment the continuous double auction of each trading day into 1440 successive 1 -second intervals. For a 
given stock, we define an average intertrade duration for each interval as follows, 



1 »n 

i = wX Tk > (3) 



where t v - is the average duration of the y'-th interval in the z'-th trading day, Nij represents the number of trades of the 
y'-th interval in the z'-th trading day. The average intertrade duration in the j-th time interval is calculated as follows, 

< T> ' = £§ T » (4) 

where Nd is the number of trading days. We note that the shape of the intraday pattern of each equity does not change 
no matter whether the durations equal to zero are considered or not. 

The results are illustrated in Fig. [2] It is interesting to observe that the two intraday patterns are qualitatively the 
same. The duration is small right after the opening of the two markets in the morning (9:30 a.m.) and in the afternoon 
(13:00 p.m.) and then increases fast to a large value followed by a relaxation in 10-20 seconds. When the time 
approaches 11:30 a.m. and 15:00 p.m., the duration sharply decreases to a small value in a few seconds. Excluding 
these short time intervals, the duration increases during the morning trading. In contrast, the duration decreases in the 
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afternoon. These patterns are consistent with those for the SZSE stocks in 2003 l26ll . The higher resolution in Fig. [2] 
allows us to uncover the fine structure in the intraday patterns. 




Figure 2: (Color online) Intraday patterns of the intertrade durations in the continuous double auction for the Bao Steel stock (a) and its warrant (b) 
traded on the Shanghai Stock Exchange in 2005-2006. 



4. Long-term correlations 



There are numerous methods proposed for the detection of long-term correlations in time series (3 
instance, researchers frequently applies the rescaled range analysis l34l - l39ll . the fluctuation analysis I4C 



For 
the de- 



trended fluctuation analysis (DFA) 141144311 . the wavelet transform module maxima (WTMM) approach 11441 14511 . and 
the detrending moving average (DMA) algorithm JHQ, to name a few. In order to investigate the possible long- 
term correlations in the intertrade durations of the Bao Steel stock and its warrant, we adopt the DFA and the DMA 
algorithms for comparison. 

4.1. DFA and DMA algorithms 

The idea of DFA was invented originally to investigate the long-range dependence in coding and noncoding DNA 
nucleotide sequences 14111 . Due to its efficiency and simplicity in implementation, the DFA is now becoming the most 
important method in the field l49ll . The DFA algorithm can be implemented as follows. 

For a given intertrade duration series {t(/)|/ = 1, 2, ■ ■ ■ , AO, define the cumulative summation series y(i) as follows, 



y(i) = Y J lT(i)-{r)], i = l,2, 



,N, 



(5) 



./=' 



where (t) is the sample mean of the r(z') series. The series y is covered by A^ disjoint boxes with the same size s. When 
the whole series y(i) cannot be completely covered by AT, boxes, we can utilize 2N S boxes to cover the series from 
both ends of the series. In each box, a polynomial trend function g(i) of the sub-series is determined. The residuals 
are calculated by 

e(i)=y(i)-g(i), (6) 
and the local detrended fluctuation function f v (s) in the v-th box is defined as the r.m.s. of the fitting residuals: 

1 vs 

[Ms)f = - X [e(0] 2 ■ (7) 

i=(v-l)i+l 
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The overall detrended fluctuation is estimated as follows 



[^)] 2 = ^-£[/vM] 2 



(8) 



As the box size s varies in the range of [20, N/4], one can determine the power law relationship between the overall 
fluctuation function F(s) and the box size s, 



F(s) 



(9) 



where H signifies the DFA scaling exponent, which is related to the power spectrum exponent 77 by 77 = 2H - 1 and to 
the autocorrelation exponent y by y = 2 - 2H ||50il5lll . 



The DMA approach is based on the moving average (MA) or mobile average technique |52ll , which was proposed 
to estimate the Hurst exponent of self-affinity signals 14611 and further developed to the detrending moving average 
algorithms by considering the second-order difference between the original signal and its moving average function 
1 47fl . Because the DMA method can also be easily implemented to estimate the correlation properties of non-stationary 



series without any assumption, it is widely applied to the analysis of real- world time series Il48ll53l459ll and synthetic 



signals as well l60ll . Extensive numerical experiments unveil that the performance of the DMA method is comparable 
to the DFA method with slightly different priorities under different situations l6ll 16211 . 

The main difference between the DFA and DMA algorithms is about the detrending procedure. In the DMA 
approach, one calculates the moving average function y{i) in a moving window J57ll . 



, r(s-i)a-fl)i 

S *=-L(i-l)9J 



(10) 



where s is the window size, |_*J is the largest integer not greater than x, \x] is the smallest integer not smaller than 
x, and 9 is the position parameter with the value varying in the range [0,1]. Hence, the moving average function 
considers f(s - 1)(1 - 6)] data points in the past and [(s - 1)9] points in the future. We consider three special cases 
in this paper. The first case 9 = refers to the backward moving average 16 ill , in which the moving average function 
y~({) is calculated over all the past s - 1 data points of the signal. The second case 9 = 0.5 corresponds to the centered 
moving average M6lll . where ~y(i) contains half past and half future information in each window. The third case 9=1 
is called the forward moving average, where ~y(i) considers the trend of s - 1 data points in the future. The residual 
sequence e{i) is obtained by removing the moving average function ~y(i) from y(i): 



e(0 = y(i) - y(f), 



(11) 



where s - [(s - 1)9] ^ i ^ N — \_{s - 1)9]. We can then calculate the overall fluctuation function F(s) as in Eq. (O 
and calculate the DMA scaling exponent through investigating the dependence between F(s) and s. 



4.2. Results 

Figure[3ja) shows respectively the dependence of the DFA and DMA fluctuations F(s) of the intertrade durations 
of the Bao Steel stock as a function of the time scale s in double logarithmic coordinates. Nice power-law scaling 
relations are observed, where the scaling ranges span over three orders of magnitude. We obtain that the DMA scaling 
exponent is Z/ s ,dma = 0.904 ± 0.008 and the DFA scaling exponent is // s ,dfa = 0.939 + 0.019, where the error bars 
are obtained at the significance level of 5%. We also perform DMA and DFA analysis on the shuffled time series of 
the intertrade durations and present the results in Fig. [3 a). We find that, for the shuffled data of the stock, the DMA 
exponent is H s , tDMA = 0.512 + 0.004 and the DFA exponent is // S ',dfa = 0.491 + 0.008. According to Fig. Oa), 
it is clear that the intertrade durations of the Bao Steel stock are long-term correlated, while the shuffled data are 
uncorrected. 

Figure[3b) shows respectively the dependence of the DFA and DMA fluctuations F(s) of the intertrade durations 
of the Bao Steel warrant as a function of the time scale s in double logarithmic coordinates. We also observe nice 
power-law scaling relations with wide scaling ranges spanning over three orders of magnitude. We obtain that the 
DMA scaling exponent is // w ,dma = 0.902 + 0.003 and the DFA scaling exponent is // Wi dfa = 0.838 + 0.017. For the 
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Figure 3: (Color online) Detection of long-term correlations in the intertrade duration time series of the Bao Steel stock (a) and its warrant (b) using 
the detrended fluctuation analysis and the detrending moving average analysis. 



shuffled data of the warrant, we find that the DMA exponent is H w > dma = 0.516 + 0.003 and the DFA exponent is 
#w',dfa = 0.503 + 0.008. According to Fig. Ob), it is clear that the intertrade durations of the Bao Steel warrant are 
long-term correlated, while the shuffled data are uncorrected. 

An interesting feature is observed in Fig.[3]that the DMA fluctuation function F(s) has a better power-law relation 
with s. Specifically, we calculate the r.m.s of the fitting errors of the eight lines and obtain that £"s,dma = 0.0397, 
£ s ,dfa = 0.0910, £ S . DMA = 0.0177, and E S ,, DFA = 0.0410 for the Bao Steel stock and £ W , DMA = 0.0143, £ w ,dfa = 
0.0785, £w',dma = 0.01 17, and £ W ',dfa = 0.0378 for the Bao Steel warrant. It is clear that 

£e,DMA < £e,DFA, (12) 

where the subscript "e" stands for s, s', w, and w', respectively. This empirical observation is consistent with the 



numerical results in Ref. [63]- However, we stress that the DMA and the DFA give qualitatively the same results. 



In the previous studies, two scaling ranges in the fluctuation function have been observed for stocks in different 



markets 11161 |24h26I1 . This crossover phenomenon disappears in Fig. [3] whose cause is unknown yet. 



5. Multifractal analysis 

Since the seminal works in the 1990's [l64U67ll . there have been a wealth of studies showing that financial mar- 
kets exhibit multifractal nature 11261 68-98]. The methods adopted in these studies include the wavelet transform 
module maxima approach Il99i - ll0lll . t he mu ltifractal detrended fluctuation analysis (MFDFA) 117311 . the partition func- 
tion and structure function approach 110211 . and so on. Here we adopt the MFDFA algorithm and the multifractal 
detrending moving average (MFDMA) algorithm proposed recently 16311 . 



5.1. MFDFA and MFDMA algorithms 

The MFDFA algorithm is a generalization of the DFA algorithm 117311 . It is necessary to note that the MFDFA 
approach was independently invented by two other groups earlier but with much less visibility in the community 
1 103 , 104 1. On the other hand, the MFDMA method is an extension of the DMA technique ||63m . which is very similar 
to the MFDFA except that they use different detrending procedures as described in Section|4] 

For both MFDFA and MFDMA, the gth-order overall detrended fluctuation in Eq. <(8j is generalized to the follow- 
ing form 



(13) 
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where q can take any real value except for q — 0. When q — 0, we have 



ln[F (j)] = — Vln[F v (j)], (14) 

■ ! „=i 

according to L' Hospital's rule. Varying the values of segment size s, we can determine the power-law relation between 
the function F q (s) and the size scale s, 

F q (s) ~ a*«>. (15) 

According to the standard multifractal formalism, the multifractal scaling exponent T(q) can be used to characterize 
the multifractal nature, which reads 

r(q) = qh(q)-D f , (16) 



where Dj is the fractal dimension of the geometric support of the multifractal measure 17311 . For time series analysis, 
we have Dj = 1. If the scaling exponent function r{q) is a nonlinear function of q, the time series is regarded to 
have multifractal nature. It is eas y to obtain the singularity strength function a{q) and the multifractal spectrum f(a) 



through the Legendre transform 110511 



a(q) = dr(q)/dq 

f(q) = qa-T(q) ' (U ' 



5.2. Results 



For the SZSE stocks in 2003, a crossover phenomenon has been reported in the detrended fluctuation functions 
lr26Tl . which is not observed in the current case for the MFDFA and MFDMA fluctuation functions. The four multi- 
fractal spectra of the intertrade duration time series of the Bao Steel stock and its warrant are illustrated in Fig. [4] It is 
found that both time series exhibit multifractal nature. 




Figure 4: (Color online) Multifractal spectra of the intertrade duration time series of the Bao Steel stock (a) and its warrant (b) obtained according 
to the detrended fluctuation analysis and the detrending moving average analysis. 



For each equity, the multifractal spectra extracted from the two methods are significantly different. The two 
multifractal spectra for the two equities obtained from the MFDMA approach are quantitatively similar, while other 
two spectra from the MFDFA method differ significantly. 



6. Conclusion 

In summary, we have performed a comparative study to investigate the statistical properties of the intertrade 
durations of a very liquid Chinese stock and its warrant using ultrahigh-frequency data. In a nutshell, we found that 
the statistical properties are qualitatively similar. The distributions of the two equities can be better described by 
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the shifted power-law form than the Weibull and their scaled distributions do not collapse into a single curve. Both 
detrended fluctuation analysis and detrending moving average analysis showed that the 1-min intertrade duration 
time series of the two equities are strongly correlated. In addition, the multifractal nature has been confirmed in the 
intertrade durations based on the MFDFA and MFDMA approaches. Contrary to the previous works, we found no 
crossover behaviors in the fluctuation functions. 
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